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We report electronic Raman scattering measurements on optimally doped YBa2Cu307_i where 
Zn or Ni impurities have been substituted to Cu. Using Raman selection rules, we have probed the 
superconducting gap in the nodal and anti-nodal regions. We show that under impurity substitu- 
tions, the energy of the anti-nodal peak detected in the superconducting state is not related to the 
critical temperature T c and that signatures of superconductivity disappear in the nodal regions. Our 
experimental findings advocate in favor of gapless arcs around the nodes. The breakdown of the rela- 
tionship between the anti-nodal gap amplitude and T c is discussed in terms of local superconducting 
gap and pseudogap. 

PACS numbers: 74.72.-h, 78. 30.-j, 74.62. Dh 



INTRODUCTION 

It is now well established that in the superconduct- 
ing state of cuprates, the gap in the excitation spectra 
has d X 2~ y 2 symmctrji This gap is usually associated 
with the superconducting order parameter A sc , and, in 
the framework of BCS theory, its amplitude is expected 
to scale with T c . An interesting way of exploring the 
link between the gap energy and T c in cuprates con- 
sists in substituting impurities on the copper sites of the 
CuC>2 planes 2 -, more specifically impurities such as Zn 
(non-magnetic) or Ni (magnetic), known for destroying 
T c without changing the doping level 3 -^. In this paper, 
we shall focus on the effects of such impurities on the 
electronic Raman response in the superconducting state 
of cuprates. 

Inelastic scattering of light by quasiparticlcs - Elec- 
tronic Raman scattering (ERS) - has the unique ability of 
probing charge dynamics in different regions of the Fermi 
surface, namely around the principal (ir, 0)(0,7r) and di- 
agonal (7r, 7r)) directions of k-space in B\ g geometry and 
B-2 g geometry respectively. In the superconducting state, 
the Big geometry probes the anti-nodal regions where 
the amplitude of the d-wave gap is maximum, and the 
B2 9 geometry probes the nodal regions, where the am- 
plitude of the gap vanishes^. Cubic and linear frequency 
dependence of the low energy Bi 9 and B2 9 responses of 
a d-wave superconductor are expected and were experi- 
mentally measured in many cuprates^i'^'^'i^. 

Previous ERS studies on magnetic impurities substitu- 
tions in cuprates (Ni in Y-123ii or Fc in Bi-2212i2) have 
shown that the Bi g peak energy seems to be insensitive 
to impurity substitutions. The case of Zn substituted Y- 
123, with only one Zn concentration, led to contradictory 
results. In one study the Bi g peak energy was found to 
follow T<~, while in another study it was argued to col- 
lapse to zeroi^ when T c = 72 K. Here, we report ERS 
measurements on a wide range of concentrations of mag- 



netic Ni and non-magnetic Zn impurities. Wc compare 
their respective effects on the Raman responses in B\ g 
and Big geometries. We show that i) the energy of the 
Bi g superconducting peak remains constant under mag- 
netic and non-magnetic substitutions up to 3%, in con- 
trast to what is expected for the pair breaking peak in 
a conventional superconducting condensate, and ii) the 
superconducting Bi g Raman response merges with the 
normal one as Zn and Ni are introduced in the Cu02 
planes. 

I. EXPERIMENT 

The synthesis protocol for YBa2(Cui_ I M 2; )307_ l 5 (M 
= Ni or Zn) single crystals have been previously detailed 
in rcf. [ll| and [la . The corresponding dc magnetization 
measurements are presented in fig. [TJ 

We have labelled the pristine crystal Y-123, the Ni sub- 
stituted crystals of ref. [ill Y-123:Ni87K (x m ~ 1%) and 
Y-123:Ni78K (xm~ 3%), and finally the Zn-substituted 
crystals of rcf. M Y-123:Zn87K (x Zn ~ 0.3%) , Y- 
123:Zn83K (x Zn ~ 0.7%) , Y-123:Zn73K (x Zn ~ 1.5%) 
and Y-123:Zn64K (xzn ~ 2%) respectively. The impu- 
rity concentrations have been checked by chemical anal- 
ysis using a Castaing electron probe, and the dT c /xz n 
(~ 15K/%) and dT c /xNi (~ 5K/%) values are consistent 
with those previously reporteds^. 

ERS have been carried out with a T64000 JY spec- 
trometer in triple subtractive configuration. Crystals 
were mounted on the copper cold finger of an He cir- 
culation cryostat and cooled down to 10 K, tempera- 
ture at which all the superconducting spectra presented 
here have been measured (100 K for the normal state 
responses). The 514 nm line of a Ar + ,Kr + laser was 
used. The laser power on the crystal surface was kept 
below 3 mW to avoid any significant heating, which was 
smaller than 3 K according to the Stokes- Anti-Stokes ra- 
tio. Big and B\ g geometries arc obtained from cross 
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FIG. 1: (color online) dc magnetization in ZFC for various im- 
purity concentrations YBa2(Cui- :E M :E )307_,5 (M = Ni (upper 
panel) or Zn (lower panel)), The Tc on set are listed in the 
table for each concentration of impurities 



polarizations of the incident and scattered electric fields 
along and at 45 degrees of the Cu-0 bounds of the CuC>2 
layers. Raw spectra I(w) have been corrected for the re- 
sponse of the spectrometer and the Bose factor n(u,T) 
to obtain the imaginary part of the Raman susceptibility 
x"(u)KJ(w)/[l + n(w,T)]. 



II. EXPERIMENTAL RESULTS AND DATA 
ANALYSIS 

A. Bi fl response of Y-123 

In figure [2]-a) are shown the Bi g Raman responses of 
the pristine Y-123 crystal in the normal (black curve) 
X^r(w) and superconducting (red curve) states. 
The normal state response consists of a flat electronic 
background with narrow phonon peaks superimposed. 
The strong peak at 340 cm -1 corresponds to the out-of- 
phase motion of 0(2,3) atoms along the c-axis. The 115 
and 145 cm -1 are respectively assigned to the motions 
of Cu and Ba along the c-axis and finally the 500 cm -1 
mode is related to the apical 0(4) in-phase motion along 
c-axis. Additional weak features close to 200 cm -1 and 
between 500 and 630 cm -1 correspond to infrared for- 
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FIG. 2: (color online)a) Bi 9 superconducting and normal Ra- 
man responses Xs( w ) ( m red)and Xjv( u) ) ( m black) of the 
pristine YBa2Cu307_^ crystal. The difference Ax"(^>) = 
Xs(w) — Xiv( w ) ( m blue) is plotted in inset, b) Raman re- 
sponses of a) after subtraction of weak phonon lines and re- 
sponses calculated using the procedure described in the text, 
c) pure electronic contributions pn and ps to the Raman re- 
sponse extracted from the procedure described in the text. 



bidden modes due to slight CuO chains disorde r 17 i 18 ' 19 . 
Entering into the superconducting state, we observe a re- 
distribution of the spectral weight from the low energy 
region to the high one (400-800cm _1 ). The weak phonon 
features are located inside the region where the back- 
ground redistribution induced by superconductivity takes 
place. According to previous studies, their temperature 
dependence exhibit small energy shifts^ (~ 2 cm -1 ) as 



well as slight broadenings (< 10%)2£, that cannot explain 
alone the renormalization effect observed here). Reso- 
nant Raman scattering measurements performed under 
ultra-violet excitation lin o 22 i 23 have recently unambigu- 
ously confirmed the electronic character of this redistri- 
bution. 

We have fitted the 100-200 cm" 1 and 500-630 cm" 1 
groups of phonon lines using simple lorentzian profiles 
in the normal state and assuming a weak temperature 
dependence, as described above, in the superconducting 
state. Figurc[2]-b) exhibits the normal and superconduct- 
ing Raman responses after subtraction of these phonon 
lines. We note that due to the 340 cm' 1 phonon, the 
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small orthorhombic distortion in and the contri- 

bution from the Cu-0 chains, the w 3 power law for the 
superconducting response of a rf-wave superconductor- 
is masked, contrary to other systems such as tetragonal 
Hg-based system o 8 ' 15 where cubic law is clearly seen. 

In order to extract the 2A value from the Raman spec- 
tra free of weak phonon lines, we use a phenomcnological 
fit of the electron-phonon coupled spectra in the normal 
and superconducting states (see refs. HH). The normal 
state susceptibility is given by : 

X;vM = Pn{uj) + + r(l + e 2 ) ^' 

where e = f2 being the phonon energy (f2~ 340 

cm -1 ), r its linewidth, and S the ratio between phononic 
and electronic matrix elements. The three last parame- 
ters are renormalized by the electron-phonon coupling V 
(a direct consequence of this coupling is the asymetric - 
Fano - lineshape of the Bi g phonon line) . 

Following the same protocol as in ref . [lj, we have accu- 
rately reproduced the electronic background in the nor- 
mal state by using the expression 

p N (u>) = C 

which is linear in frequency for tu — > and constant for 
lu — > oo (C and ldt are fitting parameters). To fit the Ra- 
man response m the superconducting state, we can 
also use eq. [T] replacing the electronic background pat(w) 
of the normal state by a renormalized pr(ui) function. 
This was achieved by adding, as in ref. ll3L two lorentzian 
profiles of opposite signs (one for the low energy loss of 
spectral weight, and one for the 400-800 cm" 1 enhance- 
ment) to pn(u>). The result of this fitting procedure is 
displayed in Fig. [2J-b). Fig. [2J-c), displays the pure elec- 
tronic contributions pjv(w) and ps(w) to the normal and 
superconducting responses deduced from these fits. 

The 2A value is assigned to the energy of the maximum 
of ps(w). It corresponds to 552 cm" 1 (8.6 ksT c ) and is 
quite similar to the energy deduced from subtracting the 
normal state contribution to the superconducting of the 
raw spectra (2A=556 cm -1 ). These values are consistent 
with those previously reported^ in YBCO system. Wc 
have also added in dotted lines on Fig. [5J-c) the calculated 
electron-phonon coupled spectra in the superconducting 
state. We notice that the energy of the maximum is 
not altered by the electron-phonon coupling and can thus 
be directly extracted from the difference x's( w ) — x'n^) 
between the normal and superconducting responses. 

B. Bi 9 response of Y-123 with impurities 

Let us focus now on the Bi g Raman spectra of Ni and 
Zn substituted Y-123 shown in Figures Q] and O respec- 
tively. The left panels of the figures show the raw Bi g 



TABLE I: Fitting parameters 



Sample 


u) (cm 1 ) 


T (cm" 1 ) 


S 


V 


Y-123 


338.6 


6.82 


-6.99 


0.79 


Y-123:Ni87K 


337.8 


9.37 


-3.86 


0.543 


Y-123:Ni78K 


339.2 


8.25 


-4.25 


0.542 


Y-123:Zn87K 


337.4 


9.27 


-2.81 


0.346 


Y-123:Zn83K 


337.1 


9.32 


-2.39 


0.330 


Y-123:Zn73K 


338.9 


9.27 


-3.28 


0.447 


Y-123:Zn64K 


338.8 


7.92 


-4.29 


0.549 



spectra, the central ones show the same spectra after 
subtraction of the small phonon lines as well as their fits 
(solid lines) with the protocol described above, and fi- 
nally, the right panels exhibit the normal pn{u) and su- 
perconducting ps(w) electronic contributions extracted 
from our analysis for each impurity concentration. 

The first striking feature is that Ni and Zn have qual- 
itatively similar effects on the Bi g Raman response of 
Y-123, with an expected^ broadening and a decrease of 
the superconducting B lg peak intensity. More surpris- 
ingly, we observe that the B lg peak energy neither col- 
lapses to zero as T c reaches 72 K (as claimed in refJ^), 
nor scales with Tj£, but remains constant over a wide 
concentration range and down to T c = 64 K. 

We have reported in figure [5] the different values of the 
pair breaking peak energies extracted from our data. The 
dashed line corresponds to the expected values of these 
energies assuming 2A/fcsX' c is kept constant. Clearly 
there is no scaling between ojb 1cj and the critical 
temperature T c . 

An independent confirmation of this fact can be in- 
ferred by looking more carefully at the fitting parame- 
ters used for the superconducting response (see table [I]). 
We find a slight increase of the Bi g phonon linewidth T 
as well as a small decrease of both the electron-phonon 
coupling parameter V and \S\ when impurities are added. 
These trends are opposite to those observed when the Bi g 
peak energy decreases and gets closer to the 340 cm -1 
phonon as, for instance, in the case of overdoping^ or un- 
der pressure^!. In these cases, spectacular changes in the 
electron-phonon coupling strength lead to an important 
increase of the phonon linewidth which is not observed 
here. 



C. B2 S response 

The nodal B 2g responses for Y-123, Y-123:Ni87K and 
Y-123:Zn83K in the normal and the superconducting 
states are displayed in Fig. [S] There are no Raman active 
phonons in this symmetry for tetragonal crystals, but due 
to the small orthorhombic distortion in the Y-123 family, 
some weak features are present in the spectra. These con- 
tributions are marginal and the main part of the spectra 
is of electronic origin. In the pristine case, we observe 
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FIG. 4: (color online) (left panel) Bi 9 superconducting and 
normal Raman responses x'sif^) ( m rec and Xjv( w ) ( m black) 
of the YBa2(Cui_ a; Ni a; )307- < 5 crystals, (central panel) re- 
sponses of the left pannel after subtraction of weak phonons 
(dotted lines) and their fits (solid lines) obtained from the 
procedure described in the text. (Right panel) Pure electronic 
contributions: ps(a>) and Pn(oj) in the superconducting (red) 
and normal state (black) respectively. 



FIG. 3: (color online) (left panel) Bi 3 superconducting and 
normal Raman responses Xs{ u ) ( m rec and Xjv(^) (in black) 
of the YBa2(Cui_ a; Zn a ;)307_5 crystals. (Central panel) re- 
sponses of the left panel after subtraction of weak phonons 
(dotted lines) and their fits (solid lines) obtained from the 
procedure described in the text. (Right panel) Pure electronic 
contributions: ps(u->) and Pat(oj) in the superconducting (red) 
and normal state (black) respectively. 



rcnormalization of the B2 ff electronic continuum in the 
superconducting state. As expected theoretically for a 
d-wave superconductor, we observe i) that the renormal- 
ization is much weaker than in the Bi s case and ii) a 
linear dependence of the B2 3 response with frequency at 
low energy. In order to extract a characteristic energy 
from our data, we have subtracted the normal state con- 
tribution from the superconducting one (see right pan- 
nel of Fig. [5]). This subtraction reveals a broad peak 
with a maximum around 470 cm -1 . It appears that the 
nodal superconducting Raman response merges quickly 
with the normal one as impurities are inserted. Indeed, 
we detect a small loss of spectral weight at low energy 
for Y-123:Ni87K (xni — 1%), but no difference between 
the superconducting and the normal i?2g responses is ob- 
served for Y-123:Zn83K (xz n = 0.7%). In contrast to 
the B lg Raman response, we cannot extract from our 
measurements a B 2 g peak energy as function of impurity 
concentrations. 
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FIG. 5: (color online) Energies of the Bi 9 superconducting 
responses (in units of T c ) under magnetic Ni (squares) and 
non-magnetic Zn (circles) impurity substitutions. 



5 



c 

3 



."3 



YB^Cb^O, 
T =92,5K 


T= 100 K 

^T = 15K 


1 1 — 

YBa,(Cu | 




T =87K(x = l%) 






T= 100 K 




T = 15 K 


1 \— 

YBa (Cu Zn) } T4 
T. = 83 K (x = 0,7%) 


1 

T= 100 K 

T = 10K 



200 400 600 800 200 400 600 800 

Raman Shift co (cm 1 ) 



FIG. 6: (color online) Superconducting and normal B2 9 Ra- 
man responses Xs( w ) ( m reci ) and Xn{u) (in black) of Y-123, 
Y-123:Ni87K and Y-123:Zn83K single crystals. The differ- 
ences Ax"(w) = x's(w) — Xn(<^) ( m blue) are plotted on the 
right panels. 



III. DISCUSSION AND COMPARISON TO 
OTHER EXPERIMENTS 

A. Nodal Response 

The disappearance of any signature of the supercon- 
ducting condensate in the the £?2g Raman response as 
impurities are inserted is consistent with previous exper- 
iments which have shown that the Big response is gener- 
ally strongly affected by any scattering source (e.g. struc- 
tural disorder) 28 i 29 . However, the effect reported here 
appears stronger than the one theoretically expected^. 
Our experimental results are in better agreement with 
the appearance of gapless arcs around the nodes in pres- 
ence of impurities for a gap of d X 2- V 2 symmetry 3 -^. In- 
deed, in this scenario, we do not expect to observe any 
difference between the superconducting and the normal 
B2 9 responses as impurities are introduced. 



B. Anti-Nodal Response 

The robustness of the Bi s peak energy when impuri- 
ties are inserted in the Cu02 plane raises the question 
whether it is a genuine result or the consequence of an 
accidental cancellation between two antagonist effects. 
Indeed a small underdoping is expected to increase the 
Bi g peak energy (ref. [3l|, H2, HI and ref. therein) while 
at the same time impurity substitutions would induce a 
decrease of the peak energy. The two combined effects 
might then result in an apparent albeit fortuitous robust- 
ness of the peak energy. 

This scenario is however highly improbable because (i) 



the crystals have been annealed under oxygen pressure 
to be optimally doped (this is confirmed by the phonon 
peak locations), and (ii) because the changes of T c due 
to impurity substitution are much more drastic in the 
underdoped side 4 -. This would imply a strong disagree- 
ment between the impurity concentrations measured by 
the Castaing electron probe and the loss of T c according 
to the dT c /x Zn (~ 15K/%) and dT c /x Nl (~ 5K/%) ra- 
tios. This is not the case here. Another possible source 
for increasing the energy of the Big peak 3 - 3 - is the pres- 
ence of an anisotropic s-wave component in the supercon- 
ducting gap. In this case a small increase of the energy 
of the superconducting peak is theoretically predicted^, 
but this effect should be rather small and is unlikely to 
cancel out the energy decrease of the A peak 3 - 4 -. 

In the absence of any other potential source for an in- 
crease of the Bi 3 peak energy we are led to conclude that 
the Big peak energy is genuinely robust under impurity 
substitutions. 

Previous angle resolved photoemission spectroscopy 
(ARPES) experiments on electron-irradiated Bi-2212 
samples (non-magnetic impurity effects in cuprates can 
be achieved by electron irradiation 3 -^) revealed a similar 
robustness of the gap energy in the anti-nodal regions. 
Consistently with our findings T c was reduced to 62 K 
and no changes in the gap energy were observed 3 -^. More 
recently, Takahashi's group has reported measurements 
of the angular dependence of the gap energy under mag- 
netic and non-magnetic impurity substitutions 3 -^ The 
gap energy in the antinodal regions was shown to re- 
main constant and gapless arcs around the nodes were 
observed, in agreement with the present findings. 

All these experiments are in agreement and they con- 
firm that the antinodal peak energy is robust under im- 
purity substitutions. This cannot be easily understood 
in the framework of conventional Abrikosov-Gorkov' the- 
ory. In this approach the presence of potential scatterers 
leads to a reduction of T c for a d-wave superconductor 
and is accompanied with a decrease of the superconduct- 
ing order parameter amplitude A 30 i 38 , leaving the ratio 
between the gap and T c constant. If we relate the energy 
of the antinodal peak to the pair breaking energy (~ A) 
as it is usually admitted^^ 8 ^^, the Bi g peak energy 
has to soften under impurity substitutions. 

At least two scenarios deserve to be explored in order 
to understand the robustness of the Big superconducting 
peak, i. e. the breakdown of the relationship between T c 
and A, under impurity substitutions in YBCO. 

The first one is built on the consensus that the Bi g 
peak is directly related to the order parameter {ojbx ~ 
2A) and based on local probe measurements showing that 
the superfluid condensate is only altered near the impu- 
rity sit o 41 i 42 , on a length scale on the order of the co- 
herence length £ (typically 15-20 A in optimally doped 
cuprates). This is the so-called "swiss cheese" picture, 
first inferred from superfluid density measurements 3 ^. As 
pointed out by Franz et al— , in this picture, the super- 
conducting gap should be considered as a local quantity 
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A(r), affected only locally by impurity scattering, and 
cannot be replaced in the gap equation by its spatial 
average, which is a crucial condition for the validity of 
Abrikosov-Gorkov' approach, ft has been shown in this 
context that the decrease of T c was faster than the de- 
crease of the average value of the superconducting gap. 
A(r) vanishes on a length scale of £ around the impurity 
site, but remains unmodified with respect to the pristine 
case far from the impurity. To our knowledge, the the- 
oretical Raman response of a superconductor with such 
inhomogeneities in the superconducting gap has not been 
calculated yet. However since the superconducting con- 
tribution to the response comes only from the regions far 
from impurity sites, we do not expect any shift of the 
pair breaking peak energy under impurity substitutions. 

The second scenario raises the question of the nature 
of the Bi g peak. It is indeed usually considered as the su- 
perconducting order parameter, but it is well-known that 
in several circumstances its energy does not scale with 
T c . This is the case in the underdoped regime where the 
B lg peak energy increases continuously as T c decreases 
(refs. I3H l32l . [33I and ref. therein) and under impurity 
substitutions, as shown by the present measurements. In 
fact, the Big peak energy behaves strikingly like the on- 
set temperature T* of the pseudo-gap. T* increases as 
the doping decreases (see for instance ref. and ref. 
therein) and is also known to be insensitive to impurity 
substitutions^^. Since near optimal doping T c and T* 
are very close, we are tempted to speculate that T* and 
u>B lg are related to the same energy scale, namely the 
pseudogap. 



Which of the two scenarios is at play in the cuprates 
is still an open question and cross-analysis between 
ARPES, Raman, tunneling and optical conductivity as 
well as further theoretical investigations are required to 
solve this important problem. 



CONCLUSION 

We have presented the superconducting Raman 
responses in impurity substituted optimally doped 
YBa2Cu307_i single crystals. For both magnetic Ni and 
non-magnetic Zn substitutions, we found that the pair 
breaking peak energy measured in the antinodal regions 
remains constant, in contrast with Abrikosov-Gorkov' 
calculations for a conventional BCS superfluid conden- 
sate. In order to explain our findings, we have explored 
two scenarios: the first one takes into account the local 
inhomogeneity of the superconducting gap and the sec- 
ond one links the antinodal peak energy to the pseudo- 
gap energy rather than the superconducting order pa- 
rameter. In the nodal regions, no signature of supercon- 
ductivity is observed in presence of impurities, suggesting 
the apparition of gaplcss arcs around the nodes. 
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